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ICCNIMUTA SPAMIIG TEEES Al® ARB0KI3SCEICES 


Abstract 

J^Lninun Spanning Trees and Minlnum Spanning Arbore- 
scence algorithns are required in mny applications, A 
nunber of algorithns are available for finding laininun 
spanning trees. This thesis presents an analytical and 
enplrical conporison of those algorithns. An algorithn 
for generating a nininun spanning arborescence, proof for 
its correctness, its conpHexity anaLysi.3, enpirical study 
and inplenentation are also presented. 



LIST OP SIIODLS 


G- Digraph 

V(g) Vertex set of G, 

L(g) Edge set of G. 

n Ctirdinality of vertex set v(g) 

e Cardinality of edge set B(g). 


A nienber of vertex set v(g) 

i^v. A nenber of edge set E(g) of a digraph G. 

X 3 

(v. ,v. ) A nenber of edge set e(g) of an undirected graph G, 

1 J 

w|(v^. }V^>| Vfeight of a directed edge^v^,v^^* 

w Kvj jV* ) I Weight of an undirected edge (v. jV. ). 

^ r j 3-3 

v/i[Gj Sim of the weights of the edges in a digraph (graph) G, 
D A Bubdigraph (subgraph) 


T A tree 


0 A cycle 

G^ The 0 reduced digraph which is sars as the givon digraph G, 

t)lbl 

G^ The i reduced digraph of digraph G^ 

G^ Big nininun incidence subdigraph cf G^* 

T- A spanning arborescence of G. 

X X 

A spanning arborescence of G. containing )E(c)j -1 edges 
fron every cycle c of G^, 



The CQ?horescence cover digraph of digraph G- 

o 

The merger tree of digrai±L 

A conposite vertex formed hy merging the vertex set of 
cycle c. 

A connected component of G^. 

Intermediate spanning’ crborescence formed by merging the 
vertices in T^ foi’ming the cycle c of G^ alone. 

Intermediate spanning arborescence formed by merging the 
vertices in forming the cycle c of G^ alone. 

Weight of the edge incident into the root of in G^, 



CHAETER 1 


IIWKODUGTIOH 

In the past four decades we have witnessed a steady development 
of grajh theoryj vAiich in the last five to ten years has blossomed 
into a new period of intense activity. Tlie main reason for this acce- 
lerated Interest in graph theory is its demonstrated applications. 

Because of their intuitive diagraiomatic representation, graphs have 
been found extremely useful in modeling systems arising in physical 
sciences, engineering and social sciences. 

"Vfaenever graph theory is applied to any practical problem, it 
almost alvays heads to large graphs - graphs that arc virtually impo- 
ssible to analyse without the help of the computer. In fact, the higii- 
speed digital computer is another reason for the recent phenomenal 
growth of interest in graph theory. 

"Although the computers are very fast, they quickly reach their 
limit if used as a brute force to solve graph theory probloms, Ihcroforo, 
an algorithm must not only solve the problem correctly, bub rausb do so 
efficiently. The two criteria for efficiency of an algorithm are the 
memory and the computation time requirements as a function of bhc size 
of the input , " 9 J • 

Perhaps the best known and most often used algorithms in graph 
theory are the spanning-tree algorithms. The spanning-tree algorithms 
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find use in various applications such as least-cost electrical wiring 
miMimum-cost conuirunication and transportation netvrorlc 
minimum cost distribution network netv/ork reliability problem. 

and the traveling salesman problem 
1 .1 WOTATlOltS iSi]D DEFIWITIOI'IS 

A directed graph or digraph is the ordered pair G- = ^v(&), e(g)> 
where 'V(g) is a finite seti|V.j,Y 2 , v^^and e(< 3') is a irreflexive 

relation in V(g). The members of v(g) are called vertices and the 
members of E, edges . The cardinality of V(g) and E{g) will be denoted 
respectively by n and e. Edge^v^jV^'^is said to be incident out of 
vertex v. and incident into vertex v.. The vertices v. and v. are 

1 ^ 1 j 

called the initial and t ermlnal vertices of edge > respectively. 

A directed graph G is said to be undirected gx'aph or simply a graph , 

if E(g) is symmetric, Edge (v. jV.) is said to be incident on vertices 

^ 0 

V. and V.* Ihe vertices v. and v. ara called the terminal vertices of 
X j X 2 

edge (v. ,v. ). 

X J 

In digre^ G, if certain members of e(g) can be placed in a sequence 

in which the first coordinate of each 
ordered pair is equal to the second coordinate of its predecessor in the 
sequence, then the set of edges in the sequence is a ycilk from v^ to v_^. 
The walk can also be written as a vertex sequence iy^, v.,v^, v^,v^). 

If in this vertex sequeiice no vertex appears more than once, the walk 
is called a path . If the initial and the terminal vertices in a walk 
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are same and no other vortex appears more than once then the walk is 
called a cycle . Vertex v^ is said to be a predecessor of vertex v^, 
if there exists a path from v. to v, and no path exists from v, to v. . 

Vertex v^ is called a successor of v. . 

^ ■ " ^ -X 

A digraph D is said to be c^'^^^graph of a digraph G- if V(d) is a. 
subset of V(g) and e(I)) is a subset of (v(d) X V(d) ) e(g) . D is said 
to be a spanning subdigraph of G if V(d) = V(g) aid is said to be an 
induced subdigraph if E(d) = (v(l)) X ^^(d) )/^ e(g) . Two subdigraphs T>^ 

Eg of a digraph G are said to be vertex disjoint if V(e^ ) r\ ^(Eg) = 0, 

®ic underlying undirected graph of a digraph G is defined as an 
ordered pear U = <V(u), E(v))> where V(u) = v(g) and (v^,v^ ) 6 e(u) if 

or i(^Vj jV^fe E(g). An undirected graph is said to be connected 
if there exists a path between every pair of vertices, A directed graph 
is said to be connected if the underlying undirected graph is connected* 

A maximl, connected subgraph of a graph is called a component of the 
graph . 

A spanning tree of graph G is a connected spanning subgraph without 
any cycles. A spanning arbor e sc enoe of a digraph is a connected spanning 
subdigraph in which every vertex other than the root has exactly one 
edge incident into it!. 

A digraph G is said to be wei/dited digrcuph, if with every edge 

<n’b>’ 

{w[<Vi,V3 


there is associated a nonnegative real number .’b>J 


for undirected graphs). 



A spanning tree T of a weighted graph is said to he a minimijm 
spanning tree (HST), if the sun of the ’./eights of the edges in the 
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tree is nininum over all the spanning trees of 

A rooted spanning tree of a y/eighted gra.ph is .said to he minimun- 
path spanning tree (mPT), if the weight of the path from root to every 
other vertex is rjininum over all such paths in the grajii. 

A spanning arhorescene of a weighted digraph is said to be a 
mininun spanning arhorescene (i®a) if the sun, of the weights of the 
edges in the arborescenco is mininm over all spanning arhorescences 
of G-. 

A spanning arborescence of .a weighted digraph is said to he nini- 
mum-path spanniiig arhorescence (iIPA) if the weight of the path from 
root to every other vertex is mininun over all diroebed paths in G, 

1 .2 SOLE BASIC PROPERTIES OP SPAMIHG IBEE 

Sone basic properties of spanning trees used in the present thesis 
are given belov/. For the proof of these properties we refer to Deo ^5 
and Harary §5^ 

theorem 1 .1 « A graph G- is called a spanning tree, if 
(i) (t IS comected and cycle-less, or 
. (ii) G is connected and has n-1 edges, or 

(iii) G is cycleless and has n-1 edges, or 

(iv) There is exactly one path between every 
pair of vertices in G, or 

(v) G is ninimally connected. 
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Kieorem, 1 ,2 > Every connected graph, has atleast one spanning tree. 

Theorem 1 .3 . In an arhoresceace, there is exactly one directed 
path from root to every other vertex. 

Theorem 1 .4 . Every strongly connected digoraph has o.tloast one 
spanning arboroscence . 

1 .3 PLAN OP THE THESIS 

Various algorithms for generation of mininua spanning trees, 
minimum-path spanning trees, minimum spanning arborescences and 
minimum-path spanning arborescence are described in Chapter 2, along 
with the modifications suggested for then, and fcheir computational 
aspects are considered. The empirical result on the computation time 
and relative merits and demerits of the algorithms are also discussed# 

On surveying the available JIST algorithms it ms observed that no 
existing I4ST algorithm could be suitably modified to yield minimum 
spanning arborescence. An algorithm for the USA is developed in 
Chapter 3. The proof of eorrfectness and computational aspects of the 
proposed algorithm are also included. 

Some applications of spanning tree algorithms are briefly described 
in Chapter 4, and some areas of future work ai'e suggested. POPiTEidT 
implementation for a few selected algorithms are given in the Appendix. 



CHAPTER 2 


smimm tree aigorithi, 2 S - a surwy 


The first ineijor contribution to minimuca spanning tree (MST) 
problem was by J.B. Kruskal [_3^* Kruskal, however, did not suggest 
any scheme for implementing his algorithm on a computer. Different 
implementations for the Kruskal algorithm have been suggested by 
0bruca|^3^, Seppanenj^4^ and Ifcllroy Another MST algorithm 

was suggested indep ardently by R.C, Prim [^2^ and E.V/. Di jkstra |j 2^. 

A recent analysis of MST algorithms by Kershunbaum and Van Slyke 
and of the set merging algorithms by Hopcroft and dllman |^2^ show 
that the Kruskal algorithm is more efficient than the Prim and 
Dijkstra algorithm vdien the input graph is sparse,, A modification to 
Kruskal 's algorithm using heap sort for picking the minimum weight 

edge and a similar modification to the Prim and Dijkstra algorithm 
using tree sort \ 3 '^ suggested by ICershanbaum and Van Siykej^JO^ 
improve the performance of the algorithms considerably. 

2.1 THE KRUSItAL AlGORITIffl 


Kruskal showed that an MST could be obtained by repeated appli- 
cations of the following steps. 

Take the smallest edge Tdiich has not been chosen or difloahdad. 
If it does not form a cycle vath some of the previously chosen edges 
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add it to the chosen edges; otherwise, discard it. 

The Kruskal algorithm is conceptually srmple, but the means of 
implementing it on a computer is not obvious. The computer imple- 
mentation vvas first suggested by A, Obruca M. In Obruca ' s imple- 
mentation the graph is represented by a nxn matrix, called the 
weight matrix , v/here the (i,j) element of the matrix represents 
the weight of the edge joining the vertices v. and v.. 

1 3 

During an iteration of the Kruslcal algorithm the subgraph formed 
by the selected edges is a collection of subtrees. In Obruca 's imple- 
mentation the selection or rejection of an edge is done using a label- 
ing procedure suggested by Doberman and Wein.b'^rger • The labeling 
procedure is as follows: 

Initialization: Initialize the predecessor of each vertex = 0. 
This makes each vertex a subtree, Dor the selection or rejection of 
picked edge (v. >v.) do the following: 

X 0 

Step 1 : If the root of subtree containing v^ the root of subtree 

containing v., then go to Step 2, else go to Step 4. 

D 

Step 2; Change the direction of the edges in the path from v^ to the 

root of its subtree. This makes v. the root of the subtree 

1 

containing it. 

Step 3- Make the predecessor of (Graft the two subtrees.) 

Step 4: Ifeke w (v. jV.) =o^r>* 

X ^ 
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Searching for the smelliest entry in the lower triangle of the * 
matrix requires computation of o(n ), The selection or rejection of 
an edge requires computation of 0(n). In the worst case the above 
steps are repeated for each element of the lov/er triangle of the 
matrix. Therefore, the execution time of Obruca's implementation is 
bounded by o(n'^). 

We now consider an implementation suggested by Seppanen [462^. 

In Seppanen' s implementation the grajii is represented by a list of 
edges (v. ,T. ), sorted in nondecreasing order of weights. During an 
iteration each subtree (i,e., connected component in the subgraph 
formed by selected edges )y. is identified by a component label. The 
labeling procedure for the four possible conditions arising in selec- 
tion or rejection of the picked edge is described next. 

Case_1_, If neither of the two terminal vertices aie included 
in a tree, the edge is taken as a new tree and vertices are labeled 
hy an increased component number C, 

Case 2 . If one terminal yertex is in a tree and the other is not 
in any tree, the edge is added to the tree. 

Case 3 . If the tvw terminal vertices are on different trees, 
then these two trees are grafted into a single tree by relabeling 


the labels of the second tree 
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Case 4 . If both the terminal v ertices are in the sairs tree, 
then this edge forms a cycle and is discarded. 

A structiJred language implementation of the algorithm is 
given next. 

Procedure SPTEEE(v,g,P,H,p,t) 

value : v,c; Inte^r v,e,p; Integer array P,H,T. 

» 

begin ! spanning tree 
c 0 
n «” 0 

for k = 1 step 1 until v ^ v(k) = 0 
for k = 1 step 1 until e do 
begin ! loop 1 
i »- P(k) 
j jS-Hffc) 

if v(i) = 0 then do 
be/^n ! 

T (k-n) fc 

if Vy ) = 0 then do 
begin ! case I; 
c c+1 
V(i) ^ 0 
y(3)'^o 
end ; case 1 ; 

else ! ! Oase 2; ); 

end ; 

if v(i) ^ v(j) then do 
bo gin ! Case 5; 
p ( k-n) k 

i V(i ) 

3 -t-- V(j ) 

for r = 1 step 1 until v do 
begin ! 

^ V(r) =r 3 then 
end ; 

end; Case 5? 
else; nie~ n+1 
end ; loop 1 ; 
p ■<- v-c+n 

end ; spanning tre e ; 
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The loop 1 of the program is executed once for every edge picked* 
Por each iteration of the outer loop one of the four cases in labeling 
procedure is carried out. Case 1,2 and 4 require a constant computa- 
tion. Case 3 requires computation linear in n, and in the worst case 
this may be executed n/2 times. Tlierefore, the upper bound on the com- 
putation is o(max(e lo^n^)). The sorting of edges requare a computa- 
tion time proportional toe log o. In actual computation it is obser- 
ved that the e log e term dominoies the computation, as the number of 
times case 3 is executed is much less than n/2. 

Ifershanbaum and Van Slyke observed that in an average case 

the number of edges picked in the I5ruslu2l algorithm is much smaller 
than e. Based on this observation they suggested a modification to 
Seppanen's implementation using heap sort . This reduces the com- 
putation time required for sorting the edges to o(e + m log e) where 
m is the number of edges picked. 

The problem of generating a spanning tree can also be formulated 
as a sot merging problem [ 2 ^ . An algorithm for grafting two sub- 
trees into a tree using set merging was suggested by Khuth ^5^ • 

The analysis of the algorithm was given by Tischer [l sj > who 
showed that o(e log e) is a lower bound, and Beterson who showed 

it to be the upper bound also. An improvement of the algorithm was 
suggested by IMcIlroy [37| . An implementation of McHroy's algorithm 


is given below. 
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Procedure SPTEEE(PjH,T,e,Ii)dge,c, w) 

Values v,c; Integers: v,e,c,Vj^jV^ lately 
Integer Array: P(c), H(c)j Edge (e), '."(e) j PredCv)} Hum(v) 
Global variables: Pred, Hun, Y,e. 



SPIEEE 

Initialization 
step 1 until v do 


begin : 
Pred(l ) ■=f~0 
Hum(l) 



ne •#«- 0 

end : Initialization 
begin : set merging 
K n e < e and c 1 then ^ 


no ^ ne+1 

0; Procedure SORT returns the subscript of the smallest v/oiglit 
edge not yot picked. 

Call SORT (W,E,Jt) 

^ p(-fc) 
y . H ( Ic) 

libeli -■f- Pind (y. ) 


labelj'^Pind (v . j 
If label i ^ labbl j then do 
begin : 

Call I/EEGE (label!, labelj ) 
end; 


end} set merging 
end} SPTEEE 



Procedure PU'ID(v^) 

Integer! top, stackp, find, label; 

Integer Array: Pred(v), IJuDa(v), stack(v); 

Global variables: Pred, I'luia, e,v; 

begin : find 
label -^v. 
stacl^ 

•yfclle pred (label) (> 0 do 
begin : 

stackp «- stackp + 1 
stack ( stackp ) label 

label ■*-- Pred (label ) 
end ; 

\d iile str?.ckp > 1 do 

begi« t 

stackp •*- stackp - 1 
top stack(staokp) 

Pred (top) label 

end ; 

end; find 

Procedure MERGE ( v. jV. ) 

Integers; idmy, v. jV.jV; 

Integer Arrays: Ered(v), Bumlv); 

Global v,ariables : Pred, Rum, v,e; 
begin ; merge 

idmy ■<^Eum(v. ) + iJura(v. ) 

K 13um(v. )> Eum(v. ) tften 
begin ; ^ ^ 

Mtimlv. ) idmy 

Pred(v. ) V. 
end ; ^ ^ 

else : 
begin: 

NumCv. ) idmy 
Pred^v. ) «- V. 
end; ^ ^ 

end ; merge. 
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An example for the above algorithm is given below. 

Let the f-H array for the inpy.t graph (figure 2,1 ) be given as 
follows : 

231 1 2325€ 

11:2354566643 5564 



n = 6 e = 14- 
Pigure 2.1 


Lhon following are the shops in the course of algorithm: 
Step 1: edge picked (l,2); selected. 



(a) 
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S-fcep 2; Edge picked (l,3)j selected 



Step 3* Edge picked (4»5); selected 



Step 4 s Edge picked (2,4); selected. 

© 

(d) 

Step 5 : Edge picked (l,5); rejected 
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Step 6; Edge picked (2,6); selected 



(f) 

Pigiore 2.2 

The algorithm termimtes at Step 6. 
Odie MST formed is 



Figure 2.3 

An analysis of the above algorithm vfas done by Hopcroft and Ullman 
I 27 ] . They shov/ed that the aLgorithm is hounded o(.e Z(n)) vAiere Z(n) 
is defined for n)> 0 to be the least number K for vtloh f(1c)^ n, viiere, 

f(o) = 1 

and F(i) = F(i-1 ^ f or i ^ 1 . 
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The first five values of 1 are 1,2,8,2024 and Z(n) is a slowly 

growing function and can be approximated to a constant. Por graphs mth 
8 to 10 vertices Z has a value between 3 and 5* 

To obtain an MST the edges are picked in increasing order of wei^t 
and they a.re selected or rejected using the Ifcllory's algorithm. Sorting 
of edges using heap sort requires conputa.tion of o(e4n log e) viiere n is 
the number of edges picked. Therefore the Kruskel ilST algorithm is 
bounded by o(e log e), 

2.2 THE PRIM MD DIJKSTEA AliOORITHM 

Tho Prim and Dijkstra MST adgorithm |l'l,4^ is based on the follow- 
irjg theorem: 

Theorem 2,1 . A spanning tree T of graph 0 is an I'tST if and only 
if for every subtree S OTj there is in T an edge of smallest weight 
among all those connecting a vertex in S to a vertex in G-S 

In the algoritlim the edges of the graph are divided into three sets. 
Set 1 t The edges assigned to the tree under construction. 

Set 2 : The edges from which liie next edge is to be added to Set 1 

will be selected (edges connecting the subtree formed by the 
edges in set I, to nei^bouring vertices). 

Set 3 : The remaining edges (rejected or not yet rejected). 

The vertices are divided in two sets. 
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Set_A: The vertices connected by the edges in Set 1 . 

B; The refflaining vertices. 

The algoritlim starts with any arbitrary vertex as set A, the starting 
vertex, and all edges incident on it as set 2. Tlie I'IST is developed by 
adding one vertex to set A in each iteration of the follovdng algorithm. 

Step 1 • The shortest edge in Sot 2 is removed from the Set 2 and 
added to Set 1 . As a result one vertex is transferred 
from Set B to Set A, 

Step 2 .* Consider the edges leading from the vertex that has just 
been transferred to Set A, to the vertices that are still 
in Set B. If the edge under consideration is larger than 
the corresponding edge in Set 2, it is rejected, if it is 
shorter, it replaces the corresponding edge in Set 2, and 
later is rejected. 

a Set 2 and Set B are enpty^etop, -othermee go to Stop 1 . 

In the implementation, the vortices ore given temporary and permnont 
labels. The vertices in Set A are given permanent labels. The vertices 
in Set B are given temporary labels denoting the weights of the edges 
(in Set 2) connecting them to vertices in Set A, The input graph is re- 
presented by an nxn weight matrix, where the(i,j)fch element of the matrix 
represents the weight of the edge joining the vertices v. and v. . An 
implementation of this is given next. 



Procedure MIHPREE (D,n, s,Pred). 

Value: ii,s; Integers: p,n,s,z 

Integer Array: D(n,n) Veet(n)} Pred(n), LabelCn) 
begin : 

for 1=1 step 1 until n do 
begin : 
label (l ) ’fr" 

Vect(l) 0 
end; 

label ( s ) '-^0 
Vect (s ) 0 

i s 
P 18 “ s 

vdille P ^ 0 do 
begin : 
n f— aft 
vect(p) 1 

p 0 

for 3 = 1 step 1 until n do 
begin s 

if vect (j ) then do 

C: Phe next statement is altered for shortest spanning tree. 

2 ^-dUjd) 

if z ~t~- label y ; then do 
begin : 

label (n ) -s-z 
Pred(j) ■*- i 

if label (j ) m then do 
begin ; 

m i!^label(3 ) 

P <-3 
end ; 
end ; 
end; 
end ; 

P 

end; spaming tree. 
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Each iteration requires computation proportional to n, and there 
are exactly n iterrutions, Therefore, the computation is bounded by ■ 



A modification to the Prim and Pijfcstra algorithm using tree sort 
was suggested by IQ3i'’shenba\mi and Van'Slyke [50]. In the modified algo-- 
rittan edges in set B are sorted , during the first iteration, using tree 
sort ^ 3^ . In every iteration, thereafter, some edges are added to the 
sorted tree (of binary tree sort), and the minimum wei^t edge is deleted 
from _ it, A POIfflEAlT implementation using the binary tree sort j^33] is ■« 
given in appendix, Ihe upper bound on computation for the modified 
algorithm pan be shovn to be o(n^) ^29 • 

2.3 WLmimi PATH SPABHIMS TBEE 

Theorem 2,2 . A spanning tree T is a, minimum path spanning tree 
of graph G if and only if for every subtree there is in'T an 

edge contained in the smallest path from the root to a vertex in G-S, 
among all the edges connecting a vertex in S to a vertex in G-S, 

Prom iHLOorcm 2.2, it is apparent that a change in the labeling 
procedure in Dijkstra's algorithm ’vrould bo sufficient to yield minimum 
path spanning tree 2J , In the modified labeling procedure a vertex 
v^ in Set B is given temporary label denoting the v/eight of the smallest 
path from root to v^, passing throu^i an edge in Set 2, 'Phis change ' 
in the implenontation of Dijkstra’s i'KP algorithm can be made by rGplacirj£ 
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■the statenent following the coinnient hy 

z ■^^label(i) + l>(i,3 )» 

Above modification requires one extra, addition for every edge in 

O 

&et 2 1ST iteration, thus requiring n additions more than the Rrim and 
Dijkstra algori'bhn described in Section 2,5. This modification, however, 
does not alter the upper’ bound on the order o-f computation or the storage 
requirement of the algorithm. ^%en Dijkstra's minimum path spanning tree 
algorithm is applied to a digraph it yields a minimum path spanning arbo- 
rescence [ 12 ] . (See I'igure 2.?). 

2.4 MlhlMOM SPAMim /JfflOEBSCEFCE 

Prom the description of MSI algorithms in Sectiin 2,1, it is obvious 
that Kruskal’s algorithm cannot yield spanning arboreacence for a digraph- 
whereas the Irin and Dijkstra algorithm will produce a spanning arbore- 
sceiice if a digraih is inputted to it. Applying the Prim and Pijkstra 
algorithm to weighted digrams one expects the resulting arboresccnce 
to be minimum spanning arborescence. But it l^ot so because Theorem 
2.1 does not hold in case of directed graphs, This can also be showir 
with the help of following counter-example. 

The given v/eighted digraph G is represented in Figure 2.4* Then 
Figure 2.5 represents the spanning arborescence T.^ obtained bj’' app:}.ying 
the Prim and Dijkstra MST algorithm and Figure 2,6 represents the 
spanning arborescence T^ obtained by applying Dijkstra 's MPT algorithm* 
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The mderlined integers in ilgures 2.5 and 2.6 represent the permanent 
labels assigned to the vertices. I'igure 2.7 represents a oinmum spa- 
nning arbore sconce for the sene digraptr. 



I'igure 2,4: Digraph G, 



Figure 2.5: Spanning Arborescence , 
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3ib,guro 2.7 ? Spanning Arborescence 


Observe that the weight of spanning arborescence In Figures 2,5 
1 

of spanning arborescencu in Figure 2.7 is 
2.5 Ji:iflP%TCAL RESULTS 

(The algorithns were coded in FOREEIAII (listings tire given in the 



and 2 , 6 are w p 


and w respectively, inhereas the woi^i 


appendix), The programs were tested on IBM 7044, mth randomly generated 

f 
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graphs of varying number of vertices and edges. The average execution 
tine for each n and e was noted using a timer routine vtoich giyss time 
in unibs of l/60 of a second. The tinier routine used takes a negli- 
gi^ble time as compared to the quantum of measurement. Tlie results are 
plotted from Figure 2,8 to Figure 2.11, 

2.6 CONCLUSION 

Speed IS not tho only measure in choosing an algorithm. Other 
important considerations arc storage requirement, form of input, availa- 
bility of algoritim and ease of implementation. In the liglit of this, 
the comparison between the various algorithms arc done. 

Tile Prim and Dijkstra algorithm is superior to the Krusknl algo- 
rithm in many of these respects. The computation time for nearly com- 
plete graphs IS small compared to Rruskal and the algorithm is ea^y to 
implement. The storage requirements are quite small, specially v/hea 
the graph is nearly complete and the edge wetighb is a simple computable 
function of the end vertices. Hero the edge weight need not be stored 
at nil, but can be oonpiitod ^/hen it is needed in Step 2 of tho algo- 
rithm. With addition of troe sort the Prim and IhQkstra algorithm 
becomes much more useful for sparse grajtis. But the algorithm looses 
its simplicity. Another advantage of the Prim and IJi^kstra algorithm 
is its adaptability to minimum path spanning tree and mnimun path 
spanning arb ore scene e. 







RATH OP COmPutATI^^N V-< OCNSITY 
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The I^ruskBl*s algorithn as iapleriented by Mbilroy is the best 
algoi^ibhiii for generating mimjaum spanning tree or miniiauia spanning 
forest for sparse graphs. 'Oaen used with heap sort, it performs well 
for density upto 0,3, The Kruskal algorithm reqmres a relatively 
large amount of memory because of the input representation of the 
graph. Also the implementation of the i5?uekal algorithm is more 
complex than the Erim and Diakstra algorithm. 



CHAITER 5 


AN ALGORITHM FOR MINIMUM SFAIttllNG ARBOIiESGEUCE 


In bhe prevxouo chapter we have discussed the various algorithms 
for generating minim'uin spanning trees. V/e have also observed that these 
algorithms cannot be modified to obtain a minimum spanning arborescenoe 
(ms a) of a given digraph. An algorithm for generating an I/ISA of a given 
weighted digraph, its implementation and computational aspects are dis- 
cussed in this chapter. 


5*1 introduction 

Let be a strongly connected weighted digi’aph. Then the minimum 
incidence subdigraph (JIIS) is defined as a spanning subdigrajii, con- 
sisting of a mimmum weight edge inoident into each verb ex of the di- 
graph G , 

0 

Let Reduced digraph G^ of G^ (or G^ \ ^ digraph obtainc 

by merging the vertices belonging to the cycles of G^ into single vertices 
Each such vertex is ceilled a composite vertex . Tho process of merging 
is defined as follows? Let c be a cycle in G^, Then all edges incident 

* 

out of any vertex v € v(o) are incident out of v , All edges incident 

X c 

into any vertex v € v(c), are incident into v , Tlio v/eights of the 
edges xnoident into v are modified as follows 5 


0 


= ’'C<vh>] - ^ 

voters ^ V(c) and ^ v^,v^^^ e(o). 


.1 
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By construc-LxoA of MIS G’, 




3.2 


Kierefore the \;ei^^hts assocx'^ted mth the edges cC redaced dig^raph 
cXco non-negative and xt Is strongly connected then is also strong- 
ly connected. 

3?o find a Laninuu spanning arboresconce of G^, we construct a se- 
quence of reduced digraphs ^m+l ^ vhere is a digraph 

v/ifch a Single vertex. This process of reduction is represented by a 
Merger Tree M. The pendant vortices and the mterriale vertices of H 
correspond to the vertices of G^ and the composite vertices at different 

stages of I’eductiohj respectively. Eoch conposlte vertex v (of G ) in 

c 1 

M has as its irmediote successors, the vertices of 6^ ^ merged to fom 

the composite vertex. The vreight of the edge connecting v to its 

c 

successor, say v , in M is wRv ,v ■0, viiere -Cy ,v > <^e(c). The ini- 

tial nnd teminal vertices of the edges corresponding to<.v y in 

digrajii G^ are also stox'^od along mth the vertex v of H* 

3 

Bob an Arborescen.ee Cpyor Di<!,raTh H be defined as spanning sub- 
digraph of G where E(h) = U U e(c) for i = 1,2,..., n. The arbore- 
o j V . nf 


£• G! 

c ^ i 


scence cover digraph contcains a uinimun sixinning arborescence. 

The algorithm for generating an MSA of a digraph G^ can now be 
stated as follows: Giv on a strongly connected weighted digrcph G^, 
obtain the merger tree M, Choose the root of the MSA as that pendent 
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vertex of norger bree vAiich. is fcirthest from the root of M, I'-arfc the 
edges of M approprio bely. The uncnrlced edges of M forn the desired 

laimuiun sprnmeg nrboresconce* 

A block digran for the algorithn is given belov/. 
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3.2 MNIMUM IMOIDSI'iCE SUBDIGRAPH 

Sone properties of m m n un mciclGnce subdi^raph are discussed 

here. 


Eroporty 3.1 . ©very vertex of a rjininuo incidence subdigreph 
has one and only one prodeces&or and 
As E(g^) contains exactly one ed£,o incident into each vortex, the 
property 3.1 is lEinediate, 

Ihoorog 3.1 . Every connected component of a minimum incidence 

subdigraph G^ has one and only one directed cycle. 

Eroof. Let Gl be a connected component in G’ , As every vertex 

in G' has a predecessor, 6' must have a cyclo. Again from property 
s s 

3.1 j jB(G^)| = h'^(^g)| f therefore, G^ cannot have more than one 

cycle . 


Using the above properties of G^ an algorithm fcac generating all 
cycles of G^, in time proportional to n can be given (see Section 3.4). 
3.3 bevel OPIilBSIT OE MSA ALGOEITHM 

Before proceeding v/ibh the algorithm, lot us prove the follovnng 
Iciins. 

Lemma 3.1 . Given any spanning ai-borescenoe of G^, there 
exists a spanning arbor escence S^ of G^ with the same root as tha.t of 
such that 

(i) contains |E(c)j -1 edges from every 
cycle c of the G^, 

(li) w [sj ^ w[a!j . 
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R?oof . If does not satisfy (i) then let contain at nost 

|e(c)| -•2 edges jCron cycle c of let €• V(c) be farthest fron 

the root (of T ) in f , such th^'t the edge <v. jV > incident into y: 

in 3?^ does not belong to E(c). lot ^ ®(c)» Then cannot 

be a successor of in T^, Eor obhervase, let be the 

path fron v- to v, in 3? , 

1 Ic 1 

By the choice of ) 

Ck-1, \/ « 


therefore, f V(c). 

the sane ergunent 2 1 ^ ^ ^ until vve coae 

to the vertex Therefore, these k vertices in the path freu to 

v^ in T^ belong to the cycle c. As is an edge in q and 

v^jVg, ... Vjj. IS a path in c fron v^ to v^. 


(^1 


(<t,,V2> , , ..., < Vi,\> . <v^^T,>}o eCo). 

Therefore |e(c)]b k or !E(c)j -1 = k-1 , which contradicts the 
assunptlon. 

Remove bhe sub arborescence rooted at v^ and abtcch it to Vj^ 

throue,h the edge <Vjj. v^> . As Vj^ is not a successor of v^ in T^, 

the resulting digraph is also a spanning nrboresoence, say T^ of 

ObeorvG thrt T’ has one nore edge fron cycle c than T has, Tlie 

edge set of T» can be obtained fron tfi-t of T by deleting the edge 
n* ” 

(v^jVi') and adding r 



equatxoia 5.1* 

w [< 

therefore , 

I ” < I" 

V<v^,v^> tE(a}^) V<v^,Vj> £E(T^) 

or 

wCp^I ^ w\.^3 3.4 

By repeating the above step at nost ^(c)| -1 tines, we obxnin 
iho spanning crboroscence T* containing |e(c)| -1 edges fron cyole c 
of As the cycles of MIS G^ are all vertex disjoint, the above 

procedure can be applied independently to each cycle of G^, On doing 
so we get the spanning arborescence of G^ containing jE(c)l -1 edges 
fron every cycle o of the MI§ G^, It is iimediate fron bhe equation 5.4 
that 

Wpj $ wy 5.5 

Q * E • D • 


GiYon a spanning nrborescenco rnd the digrajia G^, tho spanning 

arborescence can be ootained by the above lonna. Next vc shall see, 

how a spanning arborescence for reduced digraph G^_^^ of G^ can be 

obtc-ined fron 3? . 

1 
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Coa&truotion 3.1 . Lot Ijo tho given spanning nrboreaaenco of 
the digraiii G. , Bie corresponding spanning arbor esce nee T . of the 

1 1 rl 

reduced digi-'aph G of G can be obtained fi’on T as follows J 
1+1 1 1 

Step 1 . Obtain S^ froii as in Looiia 3.1 . 

Step 2: Iferge tho vertices of S corresponding to each 
cycle of G^, ^ 

Lerma 3.2 . Lot "tliG given spanning arborescence of G^ and 

bo bho corresponding spainung arborcscenco of the reduced digraph 
than 

[\+ll < "W - 

¥ c « G' 

1 

where v/fp^ is the sm of the v/ci^to of the edges in cycle o and v/e^ 

IS the weight of the edge incident into the root of in G^ . 

Proof . Lot (v^ jVg, V ,v^) bo the cycle c in the MIS G^, Note 

tliat they occure as o path say v. "v: in S (obtained fron 3? 

in step 1 of construction 5«l). Let the path denoted 

by p’. Lot V bo tho couposito vortex obtained by norglng the vertices 
0 

v^jVg,.. Vj^ , then, 

Case (i) If the set of vertices merged to fora v„ contains the 

c 

root of f (say v. ),then: let S' be an internediate arborescence obtained 
from S^ by uerging only tho vortices of cycle c, then 
E(Sj^) = e(S^) - E(P) 

B(S^) U|<vj^,v^>| - E(p) U{c-Vj^,v^>} 

= E(S^) U {< VV^>} -E(0) 
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Therefore, 

As denoted hy wo^ 

w[s^] = ''''®t “* ^IpI 

Cose 2« If the set of vertices nerged ti form v. doe.s not contain 
the root of S^, then there as an edge<JVj^,v^'^ incident into in • 
Let the edge 'v ■'/']_ correspond to the edge Then 

by equation 3.1 

» [< ^i> V] “ ” I'- "'n’’'! ^ * ['^ ^ •'^ 

IJow let be an intermediate orborcscence obtained iron hy merging 
only the vertices of cycle c, then, 

eCS^) = E(S^) D|<v^,Vj>J - ^eVV)' ■ 

= E(Sj) u^4.v^.v^>j. a -l<vy} - 

Therefore, the weight of tho arborescenco S” 

w[sj^J = w[sj - w^c] + w[^CT^,v^>] + w[<Vj^,v^‘)] - * 

Prom Equation 3.7 

w[S«] = wJSj - w[c'} 3.8 

As the cycles of 1-ffS are all vertex disjoint, merging of the ver- 
tices corresponding to any cycle c of i s independent of others. 

Therefore weight of T obtained by merging all the cycles of SHS 

It 1 

G' can bo written as follows** 

1 

^ l ® l +1 ^ "^1 

Vc C-G' 

1 


3.9 
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Fron Loniia 3«1 

Therefore 

¥c GJ 
' i 


3.10 


Q »F«X). 


'fho process of constructing fron is conbinued till we 

arrive at a spanning arlaorescenoo for the reduced digrcih 

containing only one vertex. Hero note 1h';t is a single vertex 

with no edgeo and so is S . Hierefore T . a S , Mow we shall 

n+1 m+i CL+1 

aoo how a given spanning arboresconce of can be expanded to 
obtain a spanning arborescence of • 

Oonstruotion 3.2 . Given any spanning arborescence S^, of di- 
graph G^, the corresponding spanning arborescence ^ (containing 
|E(c)j -1 edges iron every cycle c of G^ is constructed as 
follows- let y-f t'^2’ ***^ ^Ic ^ cycle, say c, of G^^^ aerged bo foru 
v^ of G^, then expand 1r^-as t'ollows ! 

Step (i ) If V is the root of S , then depending on #iich v. we 
want to be the root of delete the edge ^v^ 1 ' incident into 

the vertex in cycle c. 

Step (ii) If V is not the root S , then there is an edge <v- ,v > 
C X X c 

€ e(s ) incident inbo v.. let <v, jV^ correspond to the edge <.v ,v.'^ 

-L G >jL C I 


I 

I 


f 
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of ^ j then delete the edge m cycle c. 

Step (ill) Reassign the v/ei^t of edge /v >Y to correspond to 

'■1 c 

the weight of edge • 

By applying the above steps to each composite vertex of S^ we get 
the spanning a-rboresoence S^ ^ . 

lerim 3 .5 « Bet be the given spanning arbore sconce of and 
^ be the spanning ai’bore&concc of G^ ^ obtained by expanding the 
copiposite vortices of S^, then 


= *[\1 * 2 . 


Vc G' 
1-1 


where we^ ^ is the weight of the edge incident into the root of ^ 
in Gj_^ . 

Proof. Let S' , bo an intermediate arborescence obtained by 
1-1 

expanding the vertex v of the given spanning arborescence S , lot 
v^ be formed by merging vertices contained in cycle (v^jVg, ) 

say c of G' . Then 

Case (i) { If V is the root of S then by sicp i of oonstruction 
0 X 

3.2, the odge set of v/ith v^ as its root is given by 

E{S^__, ) = E(S^) U E(c) - (<V'^1>} * 

Therefore 

+ w[c] - w [<V^1>] 
or 

- we 


1 - 1 ' 


3.11 



Oaae (ii) J IX v 13 not the root vertex of S , then! Let the 
c i’ 

edge incident on vurtex of correspond to edge 

of G Xlien by step (ii) of construction 3 <2, the edge set of S' 

IS given as 

E(S._^ ) = E(SJ U E(o) U (<y^,v,>| - {<Vj, ,v^>} - V}' 

Iherof oro 


r 




*> - 


1 

r 

P-ij 

V/ 

s 

t ij 

+ w 

0/ 

+ V/ 

. j - " 

- w 

- I*' 


But froa Equation 3.1 

'"C<''i’Vl = ” [< wl - " 

Therefore, 


wt 




• 1 


•* 

S' . 
1-1 

Vt 

H * 
“ i 

1 b W 

^0 


3.12 


Expanding every conposite verbex in the spanning arborescence 
as above, v/e get the corresponding spanning arbor esoence , 

Then by eq'uitiom5 .11 and 12 it is iroediate that 

¥ctG' , 

1-1 

Q .E .B . 

We mil see in the follovnng theoren thob given any spanning 

arborescence of G^, tho spanning arborescence obtained by 

repeatedly locrging tho vertices of T to obtain T and then expandiiTg 

o mx I 

it to obtain S^, is the nininuQ weight spanning arborescence of G^ onong 

all the arborescencoB v/tth the s<ane root ns ® . 

o 
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^oorem 3.2 , Given any spanning arborescence and arboroscence 
cover digraph H of a strongly connoctod \/eighted digraph there 
exists a spanning arborescence of G^f vn.th the sane root as such 
that 

(i) w[sj ^ v{5j 
(ii) e(s^) c b(ii). 

Ibroof » Iiet us prove the theorem using the principle of induction* 

Given cozistruot a sequence of spanning arborescences *“** 

®D.+1 construction , 5.1 ) vhero the spanning arborescence 

consists of a dinolo vertex, say and no edges. Ib is obvious 

that root of is contained in vertex edge 

set of T , (null set) is a subset of the edge sot E of the arboro- 
n+l ^ 

scence cover aubdigraph H of G^. 

Let S . = T . , then as a basis for induction hypothesis we have 
ihtI 3n*r I 

induction hypothesis, that S^, a spanning arborescence of G^, contains 

jE(c)| -1 edges from each cycly o of G’ for o = 1,2, •».,», the root of 

S, contains the root of T„, 

1 o 

w ^ w |t J and e(S^ ) C e(H ) . 

Obtain S ' fnm S as given in Construction 5 #2 by choosing the 
i-1 3. 

root of S , , such that the root of is contained in it* S contains 
|e(c)| -1 edges from each cycle c of G^_^ and by hypothesis, 
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Gontoana |E(c)j -1 edges fron every cycle c of 0^ for 3=1, 1 +I > •••? 

(Qierefore S contcans E(c)j -1 edges fron every cycle c of G* for 

J 2 

3 5= 1-1, 1, n. Again as all the edges added to in foim^ng 

are contained in S{h), E( 3 ^ ^ ) 5 e(h), \teight of the spoiining orbore- 
scence > by Leuna 3.3 is 

Yc £ G 


Yc 6 O' 
1 


Yc € G' , 
1-1 


Proceeding simlarly till v/e arrive at , v/e nave 


= "[^n+^ -Z 


^ J, Yc^Gj 3 = i-1,i, m, 

^ — 1^1 f X ^ ^ o# 

= '"[®n+l]--Z I '■“'K 1 a 

^ ^ J, Vce<5j' . D - 

j s= i— 1 , 1, •«., n» 

By repeated applicction of Lenuo 3 »?, we have 

vtTs w(f 1 + 2 ^ 2 !. wTo]] - ‘ 2 ^ ■ 

Li-3^ LnJ j %k J- -'1^1,2, ^..,111-1. 

<1 




Yc G G» 

^ 1-1 
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therefore 

!Eherefoi-e by induction hypothesis 

and E(S ) Q E(h), 
0 


3.14 


Q •}]} *1^ • 


fhoorcu 3 .3 . Given any strongly connected v/eighted digraph 

and its arhorcscenco cover digraph H, there exists a nininuci spanning 

irrhorescencc S of G , such that S(S e(h). 
o 0 o 

Proof . Iiet bo a mninun spanning arboresoence of G^, and le t 
bo the spanning arboresoence obtained fron by Hieorea 3.1 • Then 
Vii^sJ ^ ♦ ®dt os is an tISA ^ ^iL^oJ* 

hence ~ Therefore is a mninm spanning arboresoence 

of G^ oontcaned in e(h), 

Q.B.D. 


In 'Eheoren 3.2 v/e observe that T = G and the construction 

m +1 0+1 

of G , fron G is independent of tho given spanning arboresoence T , 
u+1 o 0 

Thoroforo, given a digraph G^ alone, it is poscible to construct "r«.l 

and hence a ninimun weight arboresoence with the given root. Tlien 

it is innodiote th rt tho oi ni nun spanning arboresoence can be obtained 

by expanding the reduced digraph G^^^ = for every vertex 

as the root of S , and choosing the Euninun wei^t arboresoence among then, 
o' 
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This would require us to apply Ihe expansion procedure n tines to 
di 

reduced graph G This can be avoided if the choice of the root 

EtH 

which will yield the i4SA can be nade a pnori. The choice of the 
optiual root vertex can be nade v/ith the help of the nerger tree 
The procedure for the construction of merger tree, the spanning 
arborcscenoe from the merger tree, and the choice of optimal root 
are described below* 

Construction 5.3 . let bhe merger tree be initialised to contain 


n vertices corresponding bo the vertices of G^. let v^ bo a composite 
vertex of formed by merging the vertices of bhe cycle 
say c, of G' Tlie merging of cycle c into vertex v is represented 
by introducing vertex v into the merger tree and connecting the ver- 

G 

tices to Vq with edges ,<v^,V2>, > 

where corresponds to the edge ^ facilitate the 

expansion of to obtain S^, the initial and terminal vertices of 

edges / V A in G are also stored with vertex v of M. 

^ .L"* j IL 0 X 

The expansion of G^_^^ to obtain as given hy Theorem 3.2 can be 
implemented in terms of bhe merger treo by tho marking procedure as 
follows . 


Oonsbruobion 3.4 

Step (i)! Let V be the root of the merger tree, thon mark the 
utI 

edges in path from v^^^^ ®o* 


1 
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Step (xi) s If V IS not the root of the nerger tree, and the 
0 

edge connecting v^, to its predecessor ^corresptinding to 

edge not mrked, then nark the edge in the path froEi 

to "the nerger tree. 

Applying the above procedure breadth first to every vertex of 
nerger tree wo obtain c narlang on the nerger troo» 5lie edges of 
0-^ ^correopoioding to the unnarked edges of nerger tree represent the 
spanning arbore scenes So, vath the given root, 

Iho choice of the optmal root for con be node using the follow- 
ing bheoren, 

^eoren 3.4 . Iiot S^, with root v^, be a spanning arborescence of 
Gr^ containing j-B(o 3|-1 edges fron every cycle c of for 3 = 1 , 2 , , , .,ia. 
let M be the nerger tree obtained fron Gr^ and le t 3? denote the path 
Vj^^l , v^f fron the root of the nerger tree to then 

w[SoJ = J- w IjP]» 

i,e., the weight of the spanning arborescence so obtained by narking 
the edges of II is equal to the \reiglit of the nerger tree ninus the 
weight of the edges in the path fron root of the nerger tree to 
the root of tho spanning arborescence s^« 
li’oof t Si-on Equation 3.13 we have 



¥o-«,Gj 




« -Z 

3 = 0 , 


we 


iH 


J = 0j1 , # • . 
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As has no edges therefore 

wfs T = X wM ~ E we, 


5.15 


n 


J — 0^1^..#^ H 


Prou the construction 3.3, we have 

«M=ZZ 

) vq eG-j 

d — »•*» ^ 

v/e^ denotes tho weight of the edge incident on tho root of 
G-^. Observe that the edge incident on in merger tree also hos 
the snne weight, therefore, 


21we^ = w£p] 

3 « 0,1,,.., n 

By Equations 3*15, 15 and 17 

w[gj' = w[l-l3 - w[p] , 


3.17 


Q.E.D. 


Weight of the uerger tree w M is a constant for a given digraph 
G^, fliereforo, w[sj is nininun for thot choice of root v^ for vhich 
'^Ip} ^acooLnuni. Thus the optiaaol choice for root is thrt poxuiant 
vertex of aerger tree vdlich is farthest from the root of the merger 


tree. 
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3,4 AN MSA IMHiJjMEMTATIOH 

An Iciplenentfiti on for the construction <£ nerger tree (Ehose l), 
the opLiual choice of root vertex (Bioso II ) end the mrld.ng proceduro 
(ihase III) for obtaining an MSA of given digraph & = is given below* 
MSA ALGORlim 

Integer Axray: Mat(n,n), I/fcit i(a,n), Inod(n), ]?red(n), IiobelCn), Stack(2n), 
Mtree{2n,5), Avl(2n,2) 

Integer! n, imax, nov> cmx, cnan, rain, 1nin, iroot, tduy, id, naul, 

Paul, ipaul, troot, aroot, iroot, i,o, stack p, top, 
vertex* 

value ! n 

0 : Mot represents tho digraph in input weight natrix f oiia, 

Mati contains the imtial verta cos of every edge of 
corresponding to digroph and Matj xho terminal vortices, 
Inod gives a napping between the rov/s and oolurans, and the 
vertex of G^ it represents, Pred contains the predecessor 
list representation of the MIS and iJttiroe and nergor tree, 

Avl IS a list #j.ich contains the successor list of the 
internal vertices of merger tree. 

N is the number of verTicos in digraph G , nmx the 
roaxim.'uia nunber in merger tree at any given time and nov 


the nunber of vertices in G^j 
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begin: Initialization 

Pnia4-1 

500 * 

for 1=1 step 1 until naul do 
begin : 

Avl(i,1 0 
Avl(i,2}^i 
end; 

0 : Sot initial and torninol vortices a'^trix and imtiolizo Inod and label, 
for 1 = 1 stet) 1 unti l n do 
begin : 

Inod(i) 1 
Label (i ) 0 

Por 3=1 stop 1 until n do 
begin : 

B'lat 1 (i } 3 1 

%t j(3,i}^i 
end; 
end ; 

0 : Ilnd the minimum entry in each coli^n and subtract it fron the column. 

for 1=1 step 1 until n do 
begin ; 
nin cc 

for 3 = 1 stop 1 vintil n do 
begin : 

If I/lat (3,1 min then do 
begin : 

nin If* Mat (1,1 ) 
imin *r 3 
end; 
end; 

for 3 = 1 step 1 until n do; Ifot (3,1) ■<*-'Mat (3,1) - nin 

0 : Initialize merger tree 

Mtree (1,1 ) *<- ~1 

Mtree Ifot i(liiini) 

Mtree (i, 4 ) Mat 3 (tnin, j 
Mtree (1,5) min 

0 : Porm the predecessor list for MIS G-^« 

Pred(i ) iiain 
end ; 
nov -<r- n 
nnax •<- n 

end; Initialization, 
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0 ; Phase I; construction of Merger tree M 
0 : Part I 
begin ; 

vAiile nov 1 do 

begin s Generation of cycles of MS 
for vertex = 1 step 1 unti3 n do 
begin ; block 1 
top ^ vertex? stack p 0 
if Ifable(top) 5= 0 then do 
begins block 1 
staclep ^ stockp + 1 
stack (staclqi ) ■«- top 

Os label (top) s= 0 indicates that the vertex is encountered for the 
first tiiie. V«hen liable 5?op implies thrb the vertex has appeared 
for the second time in the path i.c» a cycle is identified, 

yMIo liable (Top) 0 do 
b eg ins 

Label (top) stackp 
staclq) ir stackp + 1 
stack, (stackp) Pred(top) 
top 'i- stack (stackp) 
end? 

end} block 1 

if Label (top) ^0 then 
begin s block 2 
cmin’^ Label (top) 
emox <- Stackp~l 
irow -le stack (cmin) 

Gall MERGE (ifet, Mabi, Mht j , Ered, Label, Inod, Miiree,n,nov, 
omax, ciain, stack) 

end} block 2 
else s 

begin ; block 3 
ygiilo staclqj''’^ 0 do 
begins 

top <- stack (stackp,) 

Label (top) ‘6' -1 
stackp •’«“ stacK> -> 1 
end} 

end} block 3 
end} 
end} 
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for X = 1 step 1 until n do 
begin t 

If Inod 0 then label (i ) 0 

endj 

end} Part 1 
0 J Part 2 

Procedure IJER&E ( Rlat, Mati, Matj, Pfed, Iiabelj Mfcree, Avl| n, nov, omln, 
cnox, stack) . 

Integer Array. Mat(n,n), Mat i(n,n), lifciij) (n,n), Pred(n)^ label (n ), 
Mti''ee(n, 5 ), Avl( 2 n, 2 ), stack (u) 

Integers: nov, onox, omn, nnax, idny, oln, Inan, ^ paTl» ipufl, 

1 , 3 , staclqp, top 
Integer value : n 

begin : 

0 : I 43 rg 0 the rows 
for 1=1 step 1 until n do 
begin : block 1 
ain 4 — Oip- 
iioin 'fr-irow 

for 3 =■ amxn step 1 until cmx do 
begin : 

idmy iS- stack(i ) 
if Mat (i,idmy) <. nan then do 
begin : 

Inin ♦-idJay 
mn l%t (i , idny ) 
end } 

irow) «.-min 
,iro\ 7 ) ^lilatj (ijinln) 
ondj block 1 

0 : IVferge the columns 

for 3 « 1 step 1 imtll n do 
begin : block 2 
•«— oo 

inin irow 

for 1 = cnin step 1 until cmax 

.ksfiia- , . 

idny stack! u ) 

3 ^ Mat (idny,o ) oin then do 
be^in: 

min ^ Mat (idny, 3 ) 
inin v- idny 


end; 
Mat(i, 
MatJ (1 
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C : llodify the predeceasor liat 

If inin = 0 and Ired(j ) = idqy then Pred(o )^irow 

end ; 

Ifo.t (irowy d ) laih 

Matx (irow, 3 ) Mat i (irow, d ) 

Matj (irow)“^ Jfet^ (iro-w^j ) 
end; 

Mat (irov/,irow) 

nmx <- noax + 1 

nov ■<?- nov - cmx + coin 


0: Pom merger tree 

for 1 = coin atop 1 mtil cinax do 
begin ; 

id ^ stack (i) 
idmy Inod(id) 

Avl(pcrfJ.,1 ) ^ idmy 
Mtree (xdoy, 2 ) ■<?' nmcoc 
ipayl «- pari 
pari V avl (pari f 2 ) 
end; block 2 
Avl(ipa'rl,2) <5- 0 
Inod (irow) nmox 

C: Pxnd the mLnxnum entry in the merged oolimn 


mxn oo 
Iron irov/^ 

Por i = 1 step 1 imtxl n do 
begin; block 3 

if Inod(i) > 0 and Mat (i, irow) mxn then do 

bogjn ; 

nin <- liilat (ijirow) 
imin i 
end ; 

end ; block 3 

0 ; Soibbxict the mnirnm entry from the colunm 

for 1 » 1 step 1 nntxl n do 
begin ; block 4 / \ /■ ^ 

if Inod(l) > 0 then Mat (x,irov/-) '^^Mht (i,±row) -• nin, 
end ; block 4 
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Oj Aasign the weight and tag to the edge incident on num: in Merger troe. 

Mbree (rmm:,5) ^ Mati(im.n, irow) 

Mtree (nEiax,4 ) -^-Mat j (lain, irow) 

Mtree (nnax, 5 ) mn 

0: Assi^ predecessor to n nnv xn G* . 

PredCirow) ^ Inin ^ ^ 

label (irow) -1 
end ; 

end} Part II » 

Of iind the optinal root vertex* 

0; Phase II 

begin ; 

troot nnox 
staolq) «- 1 

stock (staolcp ) troot 
Mbree(troot,5) ^ 0 

0; lind path length fron troot bo all teminal vertices} label vertices 
of Mtroe ; 

vtaile staokp 0 do 
begin ; block 1 
iptr-^ Mtreo (troot, 1 ) 

\’\hile iptr 0 do 
begin ; 

Idny Avl(iptr, 1 ) 
stackiJ stockp + 1 
staok(staokp) <ridmy 

Mtree (idny ) Iftree (idny, 5 )+Mtree (troot , 5 ) 

iptr •«- Avl (iptr , 2 ) 
end ; 

troot stock (staokp) 
staolqj ^ staokp ~ 1 
end} block 1 

C: PLck the vertex with largest path length (label) of Mtree 

for i = 1 step 1 mbil n do 
begin ; block 2 
If Mtroe ( 1 , 5 ) > lain then do 

begin ; 

mxn Mtree ( 1 , 5 } 
aroot ^ ± 
end} 

end} block 2 
end} Phase II 



G ; Ifergor Tree mrkirig 
G i Blase III 

begin? 

moot <- aroot 
troot «- nuax 
stnolqi *- 1 

stack (sbackii) broob 

0: Mark the edges in prbh fron troot to iroot 

vihile iroot < troot do 

begin ? 

Mtrec (iroot , 5 ) ^ -1 
iroob ^ Mbree(iroot,2) 
end; 

vAiile staclcp > 0 do, 
bo fan 

troot <- stock (stackp ) 
iptr if- BIbreo(troob,1 ) 
yjailo iptr > 0 do 

C: visit the nexrb successor 
begin? 

idoy Avl (ip br , 1 ^ 

Mbree (troot , 1 ) Avl (iptr ) 2 ) 

troot ^ idny 
staclq) staclqi + 1 
stack (staclqj ) idoy 
^ Iiltree(idU 3 r, 5)^0 bhen 

bogya- , 

iroot .«-Mtree(iday,4^ 
vAiilo xroot < troot do 
begin ? 

Mtree(iroot,5 ) ■*r -1 
iroot iftree (iroot , 2 ) 

end ; 
end ; 
end { 

stockii <^staolq)-1 
end{ 

end ; phase III 


M.T. K/'KPUR 
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C { Store the vinranrked edges in 3? oM 11 arrays 
begin > g 0 

for 1 = 1 step 1 until nmx do 
begin : 

if Mtree(l ,5) ^ 0 then do 
begin : j 3 +I 

P(j ) <- MtreG(i, 2)j H(j )^~IvflrGG(i,3 ) 

endj 

ond; 

end ; 

An Bxanulc 

Hie input digraph initial nerger tree and nininun incidence 
subdigraph are shovn in Figure 3.1, 



Input; Digraph 
Figure 3*1 
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© 00 ®® 0@0 

InliiP-1 nerger tree 



Oomposito vertex 9 Ooaposxie vertex 10 Composite vertex 11 

Iliiiiiiuri Incidence Subdigraiii 

At the end of first iteration 0. , merger tree and uininun inci- 

\ t 

dence subdigraph G| are ^own in JPigure 3.2» * 



Roducod Digraph 0^ 




iferger Tree 
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J N 



Composite vertex 12 


liLninun Incidence Subdlgraiii 0^' * 

Hgure Z.Q. - 

At the end of second iteration reduced digraxh Ggj merger tree and 
oinimun incidence subdigraph are &hov,ia in lUgure 3*5» 



Reduced Digrajh Gg 



Iferger Tree 



Oonpofntc vertex 13 


ttninm Incidence Subgraph Gg 
Figure 3.<5 ^ . 

At the cxid ol* bhird iteration reduced digr'’ph and norger tree M 
ia shown in Future 3.4. 



Reduced Digr ph 
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Phr'SG 1 of oxocirtion i& completed. Next wo noko the optmsl choice 
of the root os shown in Pigure 3.5 Jill unc’crlinc-ri intogors re- 
pro sent the woifjlit of bhc poth from idle root to the corresponding 
vertices. 



I'ho rrrldng as obtained by Hinso III is Eho’/Oi in Figure 3.C • 
odgos arc u'rlced by and the path F fron root of h to root is .sliown 

by double lines. 
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Pi£ure 3*C 


Tho lalniuvin sparining nrbore sconce forned by lAnnnrked edges of 



Ib-gure 3.7 



58 


3.5 ComTliEXlTY STUDIES 

An upper bount* on tho order of conpufation required for the 
inpXencirt'j.tioii of the mmiiun, Rpcixaiiiig arboroscence tvlgorithn is 
obtained hero. 

If’on a cursory look nt the algoribim nfc 30 apparent th-t tho 
upper bound on tho conputabi-»nal conplexity of tho ^IgoritlTn is 
deterrjined by Hinse I of tho inploucntatioHj while Ihaso II and 
Iha&o ril require tiuch less conputaxion. therefore we v/ould de'^ote 
a uaoor portion of this section m analysing Iliase I. 

I'or sinpliciby of analysis j Biaso 1 has been split into two 
pn.rts, Pert 1 gonorotes cycles of the given liHS and Part 2 nerges 
verLicee of the cycles gener-ated in Part 1 , 

Ab the ninbor of vertices or edges in any fflS does not exceed 
nj block 1 and block 3 of Part 1 roqui3?e o(n) coaputrtion per itera- 
tion, Block 2 IS oxoouted once for every cycle generated In Part 1 , 
Pus tponing the analysis of IIBliGE f or the present we analyse the 
ronainiiog portions of Phase I. 

The Blase I of the algorithn is executed a variable nunber of 
tines depending on the process of norgmg. Sach iteration reprocents 
generation of fron Gr|. Therefore, the nunber of tines Pbe&e I 

10 executed is equal to the number of ve ^’Ccg in the longest path 
fron the root of uorger tree to a pend'^^it vertex. It is apparent 
that a norgor tree (out-degree of every internal vertex being two 
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or uore) with n pendf’nt vertices \all hT.ve at nost n levels (as 
shown in Pi jure Therefore, Phase 1 ’Xiy be execute cl at 

nosi n-1 binos. 



- I Cv^l '1 

- LeorJ Z 
— L-<? v/e { 

- - Level ^ 
. L'^vel 5" 


ihjurc 

Ilonce omputi'tun tiro for Part 1 of Pnasc I is bounded by o(n'}» 
As stated above block 2 (f Part 1 is executed once for each 
cycle of BUS, As every cycle of "ICS is represented in the nercer 
tree by an interna] (couposite) vertex, block 2 is executed once 
for every internal vertex of the uerder troe« The nerjer tree 
vath n pondaifc vortices can h'^vc atnost n-1 internal vertices, the 
naxinun borne attained #ion every vertex in the tree has exactly 
two successors as shown in ]?ii 5 ures 3*»9 ^ and 3'*9« 
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Th-gore ^,9 


Tlayrurore, in the worst case block 2 is executed n-1 tires. 

Njw lot us nnrlysu the oonplexity of lart II, the IdROE procedure, 
indopendently , Hero we see thei ell four blockt are executed once 
every tiiio bho J^BRGE procedure is called, jElook 1 and BLoclc 2 re- 
quire nx 1e(c) j conputation for tacb call bo the procedure. Hence 


the total conputation tine in Mock I and 2 is proper bionnl toJ 

s: s: |l(o) 1 =n . )e{o)) 


3 0 


As can be seen f!raa the construction 5f5# sumrtion 


represents the nuiiboi* of oclgos in fclio or ^,cr troc* The liaxxoun 
Jiunber of od£;es possible in o -lort^er tree is 2n~2 (see Ihgurc 5.o^ 
Therefore, tho total conputaLion tiuo for block 1 and 2 of ?>3DE(JE 
ia at nost n(2n-2). 
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Eloclc 3 and iDlock ^ each require conputation of the order of 
and arc executed once for every cnJLl to the proco<luPGf Therefore 
total conputation line le bomded by n(n-«1 )♦ 

Hence the procedure lOSTtGE is bounded by o(n^)* Thus the order 
of conplexliy of Phase I of the algorxthn is o(u^), 

Ihoee II of tiic algorlthu is divided into tv/o blocks. Plock 1 
IS a breo traversin/j ol^joritha for the ner^jor tree, The nerger tree 
any have atuost 2n~1 vertices and 2n-2 edges, therefore, the block 1 
requires conputation of o(n). Hock 2 finds the vortex vath mxinun 
label aEiong the n teroinal vertices of the nerger tree and hence re- 
quires emputation of o(n). 

Marking procedure for the aerger tree in Ihase III is logically 
divided -nio two parto, Pirst, tree traversing which requires ccci- 
putation of the order of n, and second, edge mrlang. The nergor tree 
has atnost 2n-2 edges and atnost n-1 of those edges are mrlced. Ifork- 
ing an odgo requires a constant atiount of conputation, Therefore, 
Ruaso III of the algoT’itlin requires conputation of o(n), 

2 

Thuc conploxity of the algorithn is of the order of n , Storage 

roquirenent for the algorithn can also be seen to be of the order of 
2 

n , 

3.6 ElilPXiaCAI HESULTS AITb COHCbUSIOITS 

Tosbs sinilor to bh'^t dosenbed in CJhapter 2 v©ro carried ou' 
for the above algorithn. Tlie results are tabulated in Table 3 #1 . ♦ 



Nv n 

d\ 

10 

20 

30 

40 

50 

60 

70 

0.5 

7.0 

21 .39 

44.69 

74.3 

113.4 

161 .2 

21 6,6 

0.75 

6.8 

22.0 

43.2 

75.7 

113.5 

157.4 

212.2 

0.9 

6.8 

20.7 

42.4 

77.3 

117.7 

161,2 

207.1 

1 .0 

7.2 

20.7 

45.1 

71 .9 

116.3 

161.1 

220.2 


Table ^>1 


The variation of fche conputaLion buie with the uunbor of vertices, 
for couplotu ei^ayha la plotted rn Pieuro 3.10 J'^on i’lguro 3,10* v/e 
observe thot the avora£,o coai)utntion tine required is of bhc order of 
n * .On the avoro£,o r digraph vafch 70 vertices requires about 4 
sue aids of blue and a libble ovor 15^ words of nonory (data) locations 
It can also be seen fron the above table that the conputat:! oml true 
roquirod is indupondonb of the nunber of edges. This is so bo cause 
the input representation cf the digraph is in tho nxn weight mtrix 
forn, A li&b structure for representing the m^.ut digraph could also 
be used. But bhis would make tho algorithri nuch nore oonplo^: and 
would in fact resulb in deterioration in the average porforoance. 
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The pbove alganthn chooses the opbmnl root by itself. But 
sono of tho opplicotions aay require ,a '•Jimnun weight spajoning 
orboroscence vdth a gnvon root. This could be obt’^'inod by deleting 
Hiaso II of the algonthn and initializing taroot' to the root vortex 
desired bofire the execution of Ihasc III, If mninun rooted arbore- 
sccncG wibh different root vortices arc desired they can bo obtained 
by first generating tho uox’gei: iroo (lhase I ) and repeating the 
mrld.ng procedure (llaaso III) for every root vertex (corresponding to 
vhich niniavun arborcsconce is dosired). 



CH/ilTER 4 


COKCLUSIOIT 


Of all the analysed algorittos for finding ninanuQ spanning trees, 
tile Priri-I)i;)kstrn al;;;,onthn uodified with tree sort is the best algo- 
ritlm for nearly ccssplcbo graphs, Tlie algoritim is simple and elegant, 
a’Td can be easily nodificd to obtain mninvtia path spanning tree and 


niniraun path spanning <arboresooncG. 

Of tho diffei'ont implementations of Kruslcol's algorithn discussed, 
I'tollroy’a inpleuorrLntion ^57^ using the set merging algorithm is the 
be£Jt. 'fho IfcIlory*B implementation modiliod \iath heap sort gives 

best results for density upto 0,3. Also ’Ahon the input graph is sparse 
or when the input graph i a larger than the available memory space, iCruakal 's 
algontiim is much faster th.ti tho Prin and Dijlcstra algorithm, 

4.1 AmiOATIOMS AMD gUlTURB IROBLiiS'IS - 

Some of the applications of rumnun spanning trees were mcntioiied 
in Chapter 1 .fflxdoe applications are discussed in tnis section, 

Tho moot obvious application of I<IST is in minimum connecting net- 
vvork problem [2^* The uinimum connecting network problem can be used 
for various problems ouch os the minimum eoiuiecting road netvrark j^9j, 
uiid-mum length varing j4^ > niiiuauii cost connunioation and transportation 
networks p-1 jetc. The minimum connecting road netv/orfc 'jprolfLomp con be stated 
as follows s Suppose that vre have to connect n cities v^ , 
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sot up loast cost comodity distribution system based at ciby » 


through, a network of roads. The cost G of buildu^g a direct road 

ij 

between and is given for all peirs of cities vAiere roads can be 
bin.lt. The problem is then to find the least cost nebv/ork that connects 
all n cities. It is lEunodiately evident that this network is a miniaum 
spanning tree. 

A problem sinilax to that stated above is as follows. Suppose v/e 

have a v/arehouso at one city say v.j , and \Aant to set up a oomiaodity 

distrabution system over n cities v.^ jVg, ..., v^. The cost c^^ of 

transporting unit commodity betvraon y and v Is given for all pairs of 

^ *) 

cities which are directly coiinected to each other. The problem is to 

Ob- 
viously the solubion is a shortest path tree. This well-knov/n problem 
in Operations Resoarch la known as transportotion problem H- 

Lot us now assume thab the cost incurred in traversing any edge is 
independent of the amount carried throu^ it , In this case the solution 
is a iiininun spanning tree. 

Sene problems in social science can be modeled in the form of 

woiglitod digraphs H- Consider one such hypothetical problem of 

conslrucbing optimal command structure in a group of people, 

V . Let an edgo vV ,vN mth v/eight c denote the resi,stance of v 

to obey v . Thon the problem is to find a command structure in the 
J 

group such that the orders given by an appointed leader are best 
followed. The solution to this problem is an MSA. 
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Spanning tree algox'xthna are used ineolvin^, various graph, bheory 
problens such as finding all fundamental circmts in n graih tra- 
velling salosoan probleLi^2^j connectedness in graphs etc. One 

interesting problerj in digraphs is finding the timmum equivalent sub- 
digraph of a digraph, A oininun equivalent subdigra]ii of a digraph G 
IS a aininal subdigraph having the reachability property of G. This 
probleu ic closely relcted to the spanning tree problem. It can be 
easily seen that a spanning tree of a undirected graph G has the 
reachability property of G, But in the case of strongly connected 
digraphs this probleu can be oKovm to be atleast of the conpleocity of 
Karp-Cook class ^2^. A tight upper and lower bound on the v/eight 
of a ninlnun equivalent subdigraph of a v/oiuhted digraph and an 
approxiuato solution for it can be obtained vath the help of I ISA as 


shown boiow. 

All Imown algorithms 1 j58~| for the abo've problems use search 
bochmques. The search tine for the algorithms can be considera,bly 
reduced if bight upper and lov'or bound on 1310 \/eii_,ht of the solution 
can bo found. It is easy to see that the weight of a nimxium equiva- 
lent subdigraph is bounded bclov/ by the v/eight of the MS/i., It is 
bounded from above by the sum the v/ei^^hts of "che ninmun spanning 
arborosconcG and the nininum spanning arborescence on the digraph 
obtained from G by reversing the dii.ction of all the edges in G. 
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This bovind can bo node ti^jhter by first obtaining on MSA, replacino 
by zero the v'ei^hto of the edges of the gcoph 0- th't belong to i1S/i 
and then ob Loaning the second USA. As e natter of fact an approxi- 
uabe solution to niniuun equivalent sub-digraph is given by the union 


of bhe edges of the two nininun spanning arborescences. An exact al- 
gorithu for finding a solution to the above probleu has been given by 
Moe,lQa and Thoapson ^3^ and an efficient algorithn for acyclic di- 
graphs by Doo and Krishncunoortyj^ll^ 
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DO 120 J = O'l ,"'!UM 
JOMY=STAK(J) 
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110 I M IN = IDMY 

MIM=MAT( If lOYY) 
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33.0 CONTINUE 
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IF(XMQO( 1) )32 0,3;:0,;?^0 
330 MATUt IKDW)=viAT{ UTRDW) MIN 
B20 CONTINUe ' 

C ’ 

.;C, 3 ■"AS SIGN TMg I CH-f D;.TA(3 TO TWf f-DGF IMCIDEMT INTO 

■f;' ■ ,■;: ,,nmaM-.in th' 6 .TRfg.-, ■' 

^ x'-' 

^ ,■ 2 ;; 


^ t 'Vi 

;5* 





600 IDMY=AVLUt' r-.> ) 

STAK(STAKP)=:f DMY 
STA6P = STAKP+i 

MTP^i:E(IDMY, ^MTP^-CTnMY^-l+MTP: 

IF(AVL« IPTP,? ) HV 

61,0 IPTR = AVL( I PT’-5 , ? > 

GO TO 600 
620 STAKP = 5TAKt>»-. 

IF (STAKP)670^ e50,6,-'0 
630 IPTR=STftK( STAKPl 

TF(MTRHEnPTS tl) )620. 620,660 
640 II*IPTR 

IPTR=MTRPP( IPTR, 3 . ) 

GO TO 600 
650 CONTINUE 

PICK THE VFPT&X WITH hi?: LARGfST LAB£L. IM THE HTRFE 


STAKP=1 

STAK(STAKP) ='’mAX' 
STAKP=STAKP+i 
MAX=0 . ■ 

IMAX = 1. ■ ' 

DO 67.0 irit.M:- ■"■ 

■ 1 F "{ MTR Ei ( ■ I ■* f J -^mt, \ $ T-O ?.-6 7 Q » hb 


0 


;; v'-;. ■ 









: .'ii ^'[?v Ai 






ST AKISTAKP ) ~1 D ' Y 
IF {MTSf;S( lOMT J i ? 
721 IRCDT^NiTPFn ? DMY,4) 
GO TO. TOU 
C 

C IF NCI SUCCESS [jR L^F 

C 

730 STAKP = 5TAK.P“-i 

IF tSTAKP)790f '"90,7' 
750 1 1 =5TAK( STAKP ) 

GO TO 710 . 

790 RETURN 
END 


0 - ■'’ ? . 77 : 

■ T ij ^ i T fi. t T N : F V ~ R T ^ K 

■ 0 



EE- |v|_ 



